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Abstract 

The purpose of the present paper is to work out the general for- 
mulation of the relativistic fluid dynamics with vorticity (including 
relativistic superfluid) on the total manifold with boundary. Making 
use of the Hodge decomposition, we emphasize that the generalized 
equations of motion include a term due to non-trivial topology, while 
the pure vortex term introduced by Carter and Langlois [f?j will be 
presented merely by a co-exact form. We also consider a vortex sheet 
combined of individual vortices and discuss the boundary problem. 

1. The variational procedure yielding the relativistic fluid equations of 
motion can be discussed in diverse ways JE], 0, ||. We begin with a formulation 
which may be related both to the fluids and fields. 

Consider n-dimensional manifold X without boundary. For an exterior p- 
form a, vl ,„ Vp we define the exterior derivative D according to the formula j|, [S] 

Da =(p + 1) V [e a Ul ... Up] (1) 
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If we deal also with vectors, we can determine for an arbitrary fc-vector 
a = a Vl -- MK its divergence as [|] 

Sa= (-l)^*" 1 ) *(D*a) = V Uk a" 1 -"* (2) 

where the star operator * transforms the p-forms into the dual (n— p)-vectors, 
while * makes the inverse mapping, namely 



*a = -e^-^a = e-a .a = -r e w ... (lB o w ™"» = e • o (3) 

Here a dot denotes the inner product, e^ 1 "^™ is a fully antisymmetric tensor. 
Note that e ■ e = n and || %a = (— l) 9 ^™ -9 -* a . The divergence of an 
exterior p-form is determined as j|] 5 a = (_i) np+p+1 * (£) ^a) , where 
the Hodge operator * maps the p-forms onto the dual (n — p)-forms *a = 
i uy..v p a with the indices risen by the metric tensor . 

The Lie derivative £^ of an exterior form a along the vector field £ obeys 
the Cart an formula ^ 

£ ia = D (£■ a) +£■ Da (4) 

2. First of all, let us obtain equations of the covariant fluid mechan- 
ics on a n-dimensional space-time manifold X without boundary, dX = 0. 
Following Carter and Langlois || 0], consider a fibre bundle (X,Y,g) with 
projection g : X —>■ Y which maps the total space-time with the local 
coordinates onto (n — l)-dimensional base manifold Y with local co- 
ordinates y a where a = 1, (n — 1). A 1-dimensional fiber F y under 
the point y G Y can be interpreted as a streamline passing at this point. 
If we introduce a (n — l)-form N ai ,„ a ,~ closed on the manifold Y, the 
contra- variant functor g* : Q (Y) — > Q (X) determines a closed (n — l)-form 
iV = N ul ,, M(n _ 1} = g* {N ai _ a{n l A on X [(§]. The particle number current n 
is introduced, in light of @, as a dual to N. The closure condition DN = 0, 
according to ([I]), (|2J) leads to the conservation law V e n e = 0. 

The dynamical equations are derived from the variational formula for the 
Lagrangian 

dL= l _ H^dN^ + dw, u (5) 



2 



where the variable H = */i conjugated to AT is regarded as the momentum 
vector. (We omit the terms due to the variations of metric since they lead to 
the Einstein equations in matter not discussed here). As soon as quantum 
vortices exist in superfluid, the irrotationality condition D /j, = breaks 
down, and the Lagrangian depends on the closed vorticity 2-form w = D fi . 

Consider a fibre bundle (X, Z, p) with the projection h : X — > Z of the 
total manifold onto 2-dimensional base Z with local coordinates z (A = 
1,2). For the two-form wab closed on Z the contravariant functor h* : 
Jl 2 (Z) — > Q 2 (X) determines the vorticity 2-form w = w QV closed on X. 
The two-dimensional fibres are interpreted as world sheets of the quantum 
vortices. 

Carter and Langlois |7|] present the particle number form through a Kalb- 
Ramond field N = D B, that, in view of (0), (§), is equivalent to n e = V u b gu , 
where the bi-vector b is the dual to B. However such a representation is not 
global for a manifold X with the Betti number b n_1 (X) ^ that, indeed, 
takes place when quantum vortices appear: the closed form N, according 
to the Poincare lemma j|, is exact only on a manifold homeomorphic to a 
sphere. In the general case one should write N = F + DB, or n = f + 5b, 
where (n — l)-form F belongs to the cohomology group if ( n_1 ) (X, 9ft). 

Variations of dynamical variables in Eq. ([5]) are not arbitrary, they obey 
the convective variational principle. For an infinitesimal displacement of the 
coordinates 5y a = — £ e (dy a /dx 8 ) and 5z A = — £ e (dz A /dx e ) along an 

arbitrary vector field £ possible variations of the closed forms N and w are 
given through the corresponding Lie derivatives as 

dN = = D (i ■ N) (6) 

dw = £^w = D (£ ■ w^j (7) 

Then, substituting eqs. (§), (0) in (0) and integrating over X, we obtain the 
first variation of the action integral 

dS=JdL=(D(fN) \H) + {D(£-w) I A) = 
(-1)"" 1 (Z-F\6H) + (-l) n_1 (f • DB | SH) + (-1) 3 "" 1 (£-w\ SX) 

Let us introduce a (n — 2)-form L = *A,and present it as an orthogonal 
sum of a harmonic G, exact E and co-exact form C = 5Q by the Hodge 
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decomposition ||: 

L = G + E + C (9) 

Since vector £ is arbitrary, the eq. (H), in view of (Q), (fj) and (|3|) determines 
the equation of motion 

{f + Sb+ (-l) n+1 5c\ -w = c= *C7 

or 

{/, + V CT ^ + (-l) n+1 V ff V [(J Vv]} < = ip = *Q (10) 

In accordance to the isomorphism between the harmonic forms and the co- 
homology groups H, we may seek, without loss of generality, the solution 
as a harmonic form AF = (D 5 + SD) F = . Since /„ = (*F) V and [|J 
*AF = A * F, therefore 

V e V*f v +R*f e = (11) 

where Rf, is the Ricci tensor of the manifold X. 

So, the complete problem includes the solving of the dynamical equation 
([TOD together with eq. ([TT]) . Carter and Langlois have obtained an equation 
analogous to ([H]) on a 4-dimensional total manifold with the relevant Betti 
number being equal to zero. Now we conceive that the presence of quantum 
vortices is described by 1-form if) instead of 2-vector A — or, that eq. (|i0|) 
involves only the co-exact part of (n — 2)-form @ — , while the additional 
topological term / may exist even when the vortices are absent. 

3. The world-sheet of quantum vortices correspond to a 2-dimensional 
support SUPP w of the vorticity 2-form and, hence, to a 2-dimensional base- 
manifold Z. In principle, SUPP w can be a 3-dimensional manifold - an analog 
of the so called vortex sheet (do not mix this term with a world-sheet!) 
in superfluid He-3 0. For n-dimensional total manifold one may suggest 
SUPP w to be also a p-dimensional membrane ((2 < p < n)); this time the 
base Z will be a p-dimensional manifold. Since the vortex sheet is combined 
of individual vortices it must be an orientable manifold; the Mobius band, 
for instance, cannot be a vortex sheet. In other words a vortex membrane 
must be presented as a direct product of a (p — 2)-dimensional manifold 
and a 2-dimensional vortex world sheet. Without the loss of generality we 
can take the base manifold in the form Z = S p ~ 2 x J2, where S p ~ 2 is a 
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(p — 2)-dimensional sphere. Consider a spherical bundle (Z, J2, 7r) with the 
canonical projection 7r : Z — > X) of the manifold Z with local coordinates 
z A (A = 1, ...,p — 2) onto 2-dimensional base X with local coordinates 
a u {u = 1, 2). For an exterior 2-form closed on the reciprocal image 
Q 2 — * ^ 2 (Z) determines the 2-form w^b closed on Z. On the other 
hand, we can consider a fibre bundle (Z, S p ~ 2 , n) where the vortex world 
sheet will play the role of a fibre. Thus the closed 2-form wab is found; 
the rest procedure is similar to that discussed above. 
4. Now we consider the role of boundary. 

If the manifold X has the boundary dX, we introduce the inclusion j : 
X <— > <9X and a contravariant functor j* : f2(<9X) — > fi(X) . The tangential 
component of form a is determined as a" = while a 1 - = a — a" is 
regarded as the normal component. Integration over the total manifold with 
boundary, satisfies the formula M 

(Da, 13) x = (a, 5(3) x + ( J, ^) dx (12) 



In view of formula (|I2|) the equation (§) will contain additional boundary 
terms due to the last integral in (0). If we request the equations of motion 
to be the same as on a boundary-less manifold, these additional terms must 
vanish under the appropriate boundary conditions. Let the Betti number be 
b n_1 (X) = , that is iV = DB. Since the exact and co-exact forms belong 
to the Dirichlet and Neumann field respectively 0, we immediately write 

(Dfi) ll = (13) 

^' = {V [ ,^]} ± = (14) 

that in view of (|l|) implies no additional terms to appear in Eq. (§), and, 
hence, the equation of motion ( |IUD remains valid. When b n_1 (X) ^ , the 
particle number form N is not exact, and the extraordinary boundary term 

(N^M^gx = (*n ± ,(*ii) ± ) dx = 

under the boundary condition 

iV 11 = n 1 = (15) 



similar to (O), makes no changes in eq. (110) 
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In the non-relativistic limit we substitute n v — > n(l,v s ) , woi = , 
Wij — > ROT v s , where v s is the superfluid velocity. Then, the conditions (|15|) 
and fll4] ) reduce respectively to = and (rot-u s )" = 0, implying the ab- 
sence of flow into the walls of vessel containing superfluid and the continuity 
of vortex lines (i.e. that they are pinned to the walls or form closed rings but 
not broken up inside the fluid). 
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